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_ . Abstract 

In Maxwell's classical theory of electrodynamics the fields are frequently expressed 
by potentials in order to facilitate the solution of the first order system of equations. 
This method obscures, however, that there exists an inconsistency between Faraday's 
law of induction and Maxwell's flux law. As a consequence of this internal contra- 
diction there is neither gauge invariance, nor exist unique solutions in general. It is 
also demonstrated that inhomogeneous wave equations cannot be solved by retarded 



9i 



^ . integrals. 

g : PACS numbers: 03.50.-z, 03.50.De 
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i-£h . 1 Introduction 

Maxwell's first order system of equations in vacuo specifies the divergence and the curl 
of the electromagnetic fields E and B. In principle, this allows to calculate the fields 
■ from the given sources p and j, but the equations are coupled which confronts us with 

certain complications. In particular, it is not guaranteed that a unique solution exists 
at all, as it was questioned in [1]. The two source- free equations impose a necessary 
condition on the fields, if they exist: They must be derivable from a vector and a 
Q\ ■ scalar potential in the following way [2] : 

o: 

<N- - idA 

B = VxA, E = -V<t>--— (1) 

c at 

The potentials themselves are to be determined from the inhomogeneous equations 
which represent a coupled system of second order equations depending on the sources. 

Equation ([!]) leaves the fields unchanged when a "gauge transformation" is imposed 
on the potentials: 

A^A + ViP, (2) 

c at 

where ip is an arbitrary differentiable function. As a consequence the divergence of 
the vector potential is arbitrary to the same extent as the Laplacian Aif). This fact 
has been exploited to decouple the second order system. There is, however, no proof 
in the literature that the chosen procedure is viable and leads to unique solutions. In 
principle, the potentials are defined by the inhomogeneous equations which depend on 
the divergence of the vector potential V • A. Whether this dependence cancels when 
the solutions of the second order system are substituted into flTJ) is an open, non-trivial 
question. In Sect. 2 this problem is investigated and it is found that V • A cannot be 
chosen arbitrarily. This explains why solutions for the fields in Lorenz gauge are at 
variance with those obtained in Coulomb gauge as was found in [1] and [3]. 
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The standard procedure of solving Maxwell's equations in Lorenz gauge leads to 
decoupled inhomogeneous wave equations for the potentials which are thought to be 
solved by retarded integrals. In [1], however, it was claimed that the inhomogeneous 
wave equations cannot be solved in general, since they connect sources and potentials 
at the same time, whereas in the retarded solutions the potentials and the sources are 
to be evaluated at different times. This ambiguity is again analyzed in Sect. 3 where 
it is shown for charges moving at constant velocity that the Lienard-Wiechert scalar 
potential cannot be considered as a solution of the inhomogeneous wave equation. 

It turns out then that the system of Maxwell's first order equations does not permit 
a solution in general. Only the homogeneous wave equations, which were exclusively 
considered by Maxwell in the context of his theory of light [1], are suitable to describe 
travelling electromagnetic waves which are disconnected from their sources. The prob- 
lem is deeply rooted in an inconsistency of the first order system which is usually 
concealed by the potential ansatz (TT]). Analyzing the fields inside a plate capacitor the 
ambiguity is made visible in Sect. 4. Concluding remarks in Sect. 5 terminate this 
study on the potential method. 



2 Dependence of the fields on the divergence of the vector potential 

When we substitute the potential ansatz ([1]) into the inhomogeneous Maxwell equa- 
tions we obtain the system: 

A0=-4vrp--^ (3) 
c at 

. t* 1 d 2 A 4tt -* „ \B<f> 

where the abbreviation V • A = x was used. Naturally, the gauge function ip does 
not enter into these equations, as it cancels according to ([1]) in the expressions for the 
fields. The potentials, however, will become a function of x according to (J3J) and (@|. 
One must now investigate whether x wm also cancel in (TTJ, when the solutions of ([3]) 
and ([I]) are substituted. To this end one can exploit the linearity of eqs. (3, 4) and 
split them in the following way: 

<P = <Pi + (5) 
A<fo = -47rp (6) 

= -i£ (7) 
c at 

This set of equations is entirely equivalent to fl3]). Similarly: 

A = Ax + A 2 (8) 
. 1 d 2 A 1 4tt l„d<jn 

Applying Helmholtz's theorem on the vector potential Chubykalo et al. [4] have shown 
that the set of equations <$5§ and ([HD determines the fields uniquely, when the solutions 



4>i and A\ are substituted into ([T]). In a comment by V. Onoochin and the present 
author [5] it was pointed out that Chubykalo's procedure is equivalent to choosing 
X = 0, or adopting Coulomb gauge. It follows then by insertion of ([5]) and ([8]) into ([!]) 



B = VxA 1 + VxA 2 , E 



-V01 — V02 



c dt 



c dt 



that the terms containing \ must vanish separately 

Vxl 2 = 

c dt 

in order to render the fields independent of the chosen gauge x- 

Let us check whether the solutions of (JTj) and ffTUl) for an arbitrary choice of x satisfy 
the conditions ffl~2l) and f JTB"]) . The scalar potential 02 is given by the instantaneous 
solution of the Poisson equation ([7]) under the boundary condition 2 (66) = 0: 



(11) 

(12) 
(13) 



47T C 



d 3 x' dx(x', t) 



v \x — x'\ 



Ot 



(14) 



where the integration has to be carried out over all space. The vector potential A 2 
is given by the retarded solution of the inhomogeneous wave equation (fTUl) under the 
boundary condition A 2 (66) = 0: 



An 



An 



d 3 



x 



v \x — X 



V[x(x', t' 



i<90 2 pr, t'y 

c dt' 



(15) 



t'=t—\x—x'\/c 



The function x m ay be expressed as a product making sure that xo vanishes sufficiently 
fast at infinity: 

X (x, t) = xo (x) sin u t (16) 

From (I14D follows then: 



Ld COS Ld t 

Atc c 



X 



v \x — X 



j]Xo [X 



COS CO t 020 (x) 



(17) 



and eq. ( 1151) yields: 



A, 



47T 



d 3 x' 



uj- 



V \x — X 



— sinw (t - \x-x'\/c) V xo (x') ^-020 



Substituting the gradient of f lT7|) and the time derivative of (fTSl) into (IT5|) . and eval- 
uating at ut = 7r/2 the first term becomes zero, but the remaining volume integral 
does not vanish automatically: 



1 dA 2 (x, t) 
c 



Ot 



0J 

Att c 



d 3 x' 



v \x — X 1 



X 



sinu (\x — x 



CO 



t 1 1 



(19) 



V XO (X') - ^020 (x') 



This proves that the electric field expressed by the potentials is a function of the chosen 
divergence of the vector potential which does not cancel in (pQ). 

For the magnetic field one can draw a similar conclusion by writing the inhomoge- 
neous flux equation in integral form: 

f§.j=i(fl^ + §y d s (20) 

If E depends on x, this holds also for B due to the connection in (120]) . On the other 
hand, if one takes the curl of ([TO]) , one obtains a homogeneous wave equation for 
V x A2 which has only the solution V x A2 = assuming A2 (66) = 0. This implies 
that the magnetic field does not depend on \ in agreement with (fT2j) . but in contrast 
to (120]) . In Sect. 4 this ambiguity will be investigated in order to clarify whether 
Maxwell's equations have unique solutions at all. Before, however, let us analyze an 
inhomogeneous wave equation of type (j3J) and demonstrate in the next Section that it 
cannot be solved by a retarded integral in general. 



3 Attempt to solve an inhomogeneous wave equation 

Although the Coulomb gauge x = is the natural gauge, since it follows also from 
Helmholtz's theorem applied on the vector potential [4], most textbooks make use of 
the Lorenz gauge 

which results in an inhomogeneous wave equation for the scalar Lorenz potential by 
substitution into (13]): 



A0 L (x , t) 



1 d 2 (p L (x, t) 

? op 



-47T p (x, t) 



One may also consider the Poisson equation in Coulomb gauge 

A0 C (x, t) = -i7Tp(x, t) 

and subtract it from (122]) 

1 <9 2 0L (X, t) 



A 



(x , t) - (j) C (x, t)) 



dt 2 



(22) 



(23) 



(24) 



Note that the charge density cancels, since it is taken both in the instantaneous equa- 
tion (T231 and in the wave equation (|22p at the same time t when the potentials are 
evaluated. The formal unique solution for the difference of the potentials is the inte- 
gral: 



L (x, t) - <j) C (x, t) 



3„/ 



(fx 



-1 



\x 



1 a 2 L (f',f) 



An c 2 



dt 2 



(25) 



assuming (fic (00) = 4>l (00) = 0. This equation must be satisfied when the individual 
solutions of (122]) and (123]) are substituted. In order to check on this let us consider a 



point charge e moving along the x-axis with constant velocity v. The instantaneous 
Coulomb potential resulting from (123!) is the well known expression 



C (x,t)= (26) 
y (x — xo — v t) + y 2 + z 2 

where x ist the position of the charge at t — 0. Lienard-Wiechert have calculated the 
retarded integral for this situation obtaining [1]: 

g 

<fi L (x,t)= . r = x +vt (27) 

J(x — r) v 2 /c 2 + (1 — v 2 /c 2 ) (x — r) + y 2 + z 2 

On the x-axis expressions ( 1261) and ( 1271) are identical so that the integral on the r.h.s. 
of ( 1251) must vanish there. This is, however, not the case, if one substitutes solution 
( 1271) into (1251) and carries out the integration over all space. One obtains for r = 
(see Appendix): 

2 

4> c (x,t)- 4>l (x, t) = 6 2 V (28) 

x (c z — v z ) 

This proves that the retarded integrals are not proper solutions of the inhomogeneous 
wave equations which appear not to have a solution at all. The same conclusion was 
reached in [1] and [6]. 



4 An inconsistency in determining the magnetic field 

In order to facilitate the analysis of the flux law ( 1201) let us consider an axisymmetric 
case where a plate capacitor is charged up by a variable current (Fig. 1). In cylindrical 
coordinates one has for the Z- component of (12"U1) : 

1^ = ^+1^ (29) 
R 3R c c at 

In the region between the plates, where the conduction current vanishes, one may 
integrate ( |29l) and obtain for the circular magnetic field component 

1 f dEy 
o 

The quasi-static electric gradient field, which is created by the surface charges on 
the capacitor plates, produces obviously a magnetic field. In fact, a measurement 
of this field created by the displacement current was reported by Bartlett et al. [7] 
in agreement with Stokes law (130]) . The slope of this magnetic field was constant 
according to the results in Fig. 4 of Ref . [7] . Taking the curl of this field one obtains 
a spatially constant displacement current between the plates which is proportional to 
the time derivative of Ez- On the other hand, Faraday's law predicts for the same 
field component the expression 

ldB^ = dE 1 _dEn , . 

c dt dR dZ [ ' 



which yields upon integration over time a magnetic field exclusively created by the 
solenoidal part of the electric field. Since the irrotational quasi-static electric field 
between the plates created by the surface charges is constant in space as measured 
in [7], the r.h.s. of ( 13~Tj) vanishes and it would be impossible to build up a magnetic 
field between the plates. Nevertheless, in agreement with (130 p the changing irrotational 
gradient field is well capable of producing a magnetic field in contrast to the prediction 

of dSU). 




Figure 1: Magnetic field created by a quasi-stationary electric field in a plate capacitor 

There is apparently an intrinsic inconsistency between Faraday's law of induction 
and Maxwell's flux law. The discrepancy is not obvious as long as the potential ansatz 
([I]) is adopted. It guarantees that Faraday's law is satisfied automatically once the 
vector potential is determined from the flux law. A temporal evolution of the magnetic 
field, however, is only possible in the presence of a solenoidal electric field E s according 
to Faraday's law: 

B = -c [ dtV x E s (32) 



whereas by spatial integration one finds that the magnetic field is also a function of 
the changing irrotational electric field E t produced by charge separation: 

In general, equations ( 1321) and ( 1331) are incompatible. 




5 Concluding remarks 

The analysis presented in this paper forces us to recognize that Maxwell's system of 
first order equations cannot be solved consistently as it contains an internal contra- 
diction. The potential method - which was already adopted by Maxwell himself - 
conceals this fact, but it allows deriving a second order wave equation for the vector 
potential modelling electromagnetic waves successfully. Maxwell considered the homo- 
geneous wave equation in a region far away from the sources and formulated boundary 



conditions for its solution [1]. Possibly, he was aware that an inhomogeneous wave 
equation is not solvable by retarded integrals as shown in Sect. 3. This is also true 
for the vector wave equation (j3J) regardless which gauge is chosen. 

Classical electrodynamics requires apparently a thorough revision with special at- 
tention to the interaction of waves with matter. At this point a concrete proposal is 
not available, but it is likely that Planck's constant must be built into the system of 
equations, since it plays a major role in the quantum theory of light that has replaced 
Maxwell's theory of light. 



Appendix 

In Sect. 3 the volume integral (125]) had to be evaluated: 



d 2 <p L (£' , t) 



d 3 x' 



An c 2 

e/3 2 



An c 2 J J J \.r — 

i o 2 



dt 2 



d 3 x' 



\x — x' I dt 2 



x' - x - v t) 2 + (1 - (5 2 ) (y' 2 + z' 2 ) 



An 



3 ' 



d 6 X 



-2 (x' - x - v tf + (1 - (3 2 ) (y' 2 + z' 2 ) 



{x> - x - v t) 2 + (1 - /? 2 ) (y> 2 + z' 2 ) 



where (3 = v/c. Changing variables one may write 

Evaluation at the time t = —x /v when the charge has reached the origin yields on 
the x-axis: 



An 



d 3 x' 



-2 (^ + ^ + (1-^) (R 2 + R: 



(x + R x f + (1 - (3 2 ) (m + Rl 



In spherical coordinates one has 

R x = R cos 6 , R y = R cos ip sin 9 , R z = R sin ip sin i 
and the volume element becomes 

d 3 x' = R 2 sin 6 dp d6 dR 

This results in: 



2tt 



RdR / sinOdO / dp 



I= e Jl 
An 

oo o 
Integration over the angles yields 



/ = - e (3 2 x J RdR 

Q 



-2 (x + R cos Of + (1 - f3 2 ) R 2 sin 2 ■ 
[x + R cos6) 2 + (1 -P 2 ) R 2 sin 6 



(R - xY (R + x)-^/(R + x) z (R-x 
(R-x) 2 J(R+x) 2 (x 2 - R 2 p 2 Y 



For R > x the integrand vanishes, and for R < x the integral becomes 



X 



xR 



-e/3 2 



I = -2e(3 2 J 



dR 




R 2 f3 2 ) 2 



x (1 -(3 2 ) 



o 
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